Observation of free hole gases at ambient conditions 
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By studying fluctuations in the electrostatic potential ol a single electrode, we were able to 
perform in situ electron spectroscopy. Electron exchange processes that occur at the surface of the 
electrode were triggered by the surface Auger process either using metastable molecular nitrogen or 
excited state helium at ambient conditions. Ongoing redox reactions were decomposed into their two 
fundamental half reactions and the reduction and oxidation potentials were detected. We found that 
redox reactions near the electrode surface are the result of binary interactions of a free electron gas 
with a free hole gas that occur on a femtosecond time scale. The measured lifetimes of electron-hole 
recombination processes and the elastic scattering rate of the hole-hole process are in fair agreement 
with theoretical estimates. The observed asymmetry in the energy distributions of free electron and 
free hole ensembles suggests that at thermal equilibrium the identity of hole differs from that of an 
electron vacancy. 

PACS numbers: 82.80.Dx, 52.25.Kn 

Keywords: electron ensembles, thermal equilibrium, statistical mechanics, quantum fluctuations, Wiener- 
Kinchine, electron spectroscopy 
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I. INTRODUCTION 

The ultimate goal in quantum chemistry is to achieve 
control over molecular reactions at a single electron 
level. The development of diagnostic tools to detect spe- 
cific electron processes has been a continuous effort for 
decades Q], [H, y] . Of practical significance are those tools 
that can operate under ambient conditions such as in situ 
STM and electrochemical devices @, H, @] . 

The physics and chemistry studied in a system that is 
held under ambient conditions (high particle densities), 
is fundamentally different from the physics and chemistry 
of the same system in a vacuum (low densities) . Under 
conditions where a free molecular flow is established, no 
scattering occurs until a reaction takes place. At high 
particle density elastic scattering rates are such Q , that 
a complete thermalisation of the system has occurred by 
collisions of electrons with molecules in the background, 
before and after reactions have taken place. As a result, 
the response of an electron system that is evaluated at 
high densities proceeds through the response of an ensem- 
ble of particles that has reached a thermal equilibrium. 
This is in contrast to the low density case where the re- 
sponse of the system is evaluated by the response of a 
single electron at a fixed energy. 

Electrochemical devices are based on the detection and 
quantification of charge transfer. In the analysis of redox 
reactions the energy required in the process of electron 
transfer is a measure of the standard redox potential of 
the specific molecule under investigation. If one is able to 
quantify electron energies with sufficient resolution it is 
possible to discriminate between different ongoing redox 



'Electronic address: hennie.mastwijk@wur.nl 



reactions and resolve the molecular chemistry to the level 
of single electronic quantum states 0, 0] . In this paper 
we provide experimental evidence of the possibility of sin- 
gle electrode quantum coulometry as tool to perform in 
situ high resolution electron spectroscopy. In essence we 
present a method which allows the study of isolated half 
reactions that corresponds to well defined reduction (oxi- 
dation) states by a direct count of the number of electron 
(hole) states and differentiate with respect to their redox 
potential. An important feature of the method is that 
measurements are carried out at ambient, nearly field free 
conditions without the principle requirement of a counter 
electrode. In order to understand the electronic response 
of a single electrode subjected to chemical processes at its 
outer surface, we have re-examined the physics of electron 
ensembles at the interface of an electrodeQ using the 
principles of surface Auger spectroscopy [!|. We found 
that neither high vacuum conditions nor electron spec- 
trometers employed in conventional Auger spectroscopy 
are requisites. The entire experimental set-up was re- 
duced to a single working electrode operated at ambient 
conditions, directly recording the redox reactions near 
the interface which were detected as events that cause 
fluctuations in the electrode potential. 

This paper is organized as follows. First, we shall dis- 
cuss principle features known in quantum physics con- 
cerning the measurement of electronic energies and moti- 
vate how voltage fluctuations are converted into electron 
energy spectra. Secondly, we shall provide experimen- 
tal evidence that electron energy spectra can be derived 
from fluctuations in the electrostatic potential by study- 
ing the surface Auger process invoked by excited state 
nitrogen and helium when these gases are exposed to a 
charged metal surface. Finally, we will argue that redox 
events are the result of electron-hole recombination in- 
teractions of a free electron gas and a free hole gas that 
are omnipresent. 
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II. THEORETICAL FRAMEWORK 

The electronic wave function of a discrete quantum 
system is expressed in terms of its eigen states tpi by 

* = Clip! + C 2 ^ 2 + -Cntpn (1) 

The expectation value of the energy is given in terms 
of its eigen energies Pj by [l(| 

(E) = c\Ei + c 2 2 E 2 + ...c 2 n E n = P,E t (2) 

where Pj is the probability of finding the system at 
a specific eigen energy Pj. It is important to realize 
that the expectation value is not a true observable but 
a derived quantity. Certainly, a detector can be config- 
ured to perform an averaging procedure to read an actual 
outcome of magnitude (P) but the only irreducible ob- 
servables (actual outcome of an measurement) are the 
eigen energies Pj. As an example, consider playing dice. 
The outcome of each throw is an element of the set 
(1,2,3,4,5,6). After collecting a large number of outcomes 
the expectation value, or arithmetic mean, approaches 
the value 3^. However, the expectation value is never 
observed as an actual outcome. 

The determination of the spectrum of eigen energies 
is considered in this paper to be the result of a stochas- 
tic process. By repetitive measurements, a probability 
Pi can be assigned to find the system at a specific eigen 
energy Pj. In an experiment this is achieved by counting 
the number of (uncorrelated) occurrences of Pj. In the 
discrete case the relative frequency of occurance (Pj) of 
eigen energies Pj determine a discrete spectrum whereas 
in the continuous case the values of Pj are replaced by a 
density of states fiol ] . For a quantum mechanical ensem- 
ble in thermal equilibrium the probability Pj is propor- 
tional to D(Ei) where D(E,kT) is the proper quantum 
distribution function of energies of the ensemble [T3| un- 
der consideration at thermal energy /3 = kT under con- 
ditions of thermal equilibrium. The energy distribution 
dn{E)/dE of the ensemble is defined by the product of 
the single particle density of states p{E) and the quantum 
distribution function 0, Ho| . For a dense, free fermion 
ensemble the energy distribution is according to: 



fermion ensembles. We will use eqn. |3]to obtain the tem- 
peratures of electron and hole gases that are contained 
near the surface of a charged metal electrode. 

The Fermi level of an ensemble of electrons is de- 
termined by the energy of the two last added particles 
[1] • Since the energy of a single electron is given by 
P = -qV [IH, where V is the static electric potential and 
q = 1.602 10~ 19 C, the Fermi energy of an electron ensem- 
ble is given by Ej = —qV. Here we make use of the fact 
that the static electrical potential of an metal object is a 
measure of the Fermi level of the electron ensemble that 
is confined at its outer surface. The quantum mechanical 
response of the Fermi energy upon exchange of a single 
electron is determined by the density of states of surface 
charge (quantum capacitance per unit of area) given by 
Cqm = q 2 m/nh 2 = 4e /a Q where a is Bohrs' radius. 
The change in the electrostatic energy of a macroscopic 
amount of charge (ensemble) is determined by the clas- 
sical density of states of surface charge (capacitance per 
unit of area). For a (macroscopic) capacitor with radius 
r this is given by c c / = e /r [111 ]. Thus, the change in 
energy by the release of a single electron charge (quan- 
tum mechanical regime) is identical to the release of a 
macroscopic amount of charge by an ensemble (classical 
regime) of magnitude Q = Ar/a a . 

To obtain an electron energy spectrum rather than an 
expectation value, we detect fluctuations in the static 
potential of the metal electrode rather than averages. 
The electron energy distribution is obtained by repetitive 
sampling of uncorrelated readings of the static electrode 
potential. Random fluctuations in the energy of a sys- 
tem at stationary conditions are assumed to obey the 
Wiener-Khinchine relations 0. The correlation func- 
tion K(s) of random temporal fluctuations in the en- 
ergy E(t) of an ensemble is given by the self convolu- 
tion K(s) = (E(t)E(t + s)) 0. This correlation function 
is the Fourier transform of the (power) spectral density 
of E[t) which is Pj for a discrete quantum system and 
dn(E)/dE for a continuous system. The spectral density 
dn(E)/dE can be readily obtained in an experiment us- 
ing the following procedure: the number of samples dn 
that is found in the interval —q(V, V + dV) are assigned 
to the energy distribution function dn(E)/dE in the in- 
terval (P, E+dE). In other words: to obtain the electron 
energy distribution we build a histogram from observed 
series of fluctuating voltages. 



dn(E) , . , . \/2to 3/2 /= 1 
-^=KW) = ^ ^ (g _ g/)+1 (3) 

where m is the mass of a single particle in the ensemble 
and Ef the Fermi energy 0, S, The latter equation 
is our model for the energy distribution of free electrons 
(dn e (E)/dE) for a systems with total charge Q < or 
the energy distribution of free holes (dn,h(E) / dE) for a 
system with Q > 0. The characteristical shape of the 
energy distribution allows to make a comparison of ex- 
perimental data within the theoretical framework of free 



III. EXPERIMENT 

The experimental set-up of the detector that is used 
in this work is depicted and explained in figure [1] A 
Cu wire (30xlmm) serving as an electrode is configured 
as a large bandwidth, high impedance voltage probe (see 
Appendix). Its electrical potential is perturbed (figure [3]) 
either by excited state molecular nitrogen or excited state 
(atomic) helium produced in a cold plasma operated 
at ambient conditions. Specific features concerning the 
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FIG. 1: Overview of the experimental set-up. The static elec- 
tric potential of a Cu electrode is perturbed by the afterglow 
of a stationary jet of excited state helium or nitrogen pro- 
duced by a high voltage discharge. The electrode is located 
25 mm above the gas outlet. A grid near the Cu electrode, 
held at a fixed potential, is used to remove space charge sur- 
rounding this electrode. The system is operated at ambient 
conditions. Events that change the potential of the electrode 
are recorded by a digital oscilloscope as changes in the voltage 
V(t). 



afterglow of nitrogen (Vegard-Kaplan system) are dis- 
cussed in more detail in the Appendix. A preset bias 
voltage is used to add electrons (Q > 0) or holes (Q < 0) 
to the ensemble held at the surface of the Cu electrode. 

Redox events at the electrode surface involving the ex- 
change of charge, lead to a temporary change in the static 
potential of the electrode. The corresponding changes 
in voltage were recorded. The energy distribution of 
either electron or hole ensembles were recovered from 
these recordings by building a histogram as discussed 
above. Voltages were typically sampled at a rate of 1 
kHz. Excessive space charge (figure [5]) that builds up 
during the experiment was removed using an alumina 
grid (20x20mm, mesh 3mm, 70% open). 

The observed events, fluctuations in the voltage, are 
triggered by the Auger process at the metal-air interface 
when excited state nitrogen molecules are de-excited by 
surface electrons upon approach [9j as depicted in fig- 
ure [21 Note that by applying this method we directly 
detect the energy distribution of the electron ensemble 
near the electrode rather than the energy state of a sin- 
gle electron that is first extracted and then analysed on 
the basis of its kinetic energy (under vacuum conditions) 
S El- This feature enables operation at ambient gas 
conditions. At a positive biased detector a similar pro- 



cess occurs that results in the detection of the energy 
distribution of positive charges (holes). 

The assumed relation between the single particle elec- 
tron energy and the changes in the static potential of 
the detector, as discussed in section UH requires exper- 
imental proofing. We have validated that the observed 
random fluctuations in the static potential directly cor- 
respond to the electron density of states using the He I 
system. This was achieved by replacing the carrier gas 
N2 by He (Hoekloos-Linde, gas purity grade 5.0). The 
primary content of the helium plasma that was obtained 
is singly excited, triplet helium as indicated by the light 
emission spectrum given in figure SJA. This was confirmed 
by the measurement of the electron density of states de- 
picted in figure [3J3. The electron energy distribution was 
obtained at a bias of -50V by analysis of the fluctuations 
in the electrostatic potential of the Cu electrode similar 
as observed in the case of nitrogen and similar as dis- 
cussed in section |Hj The electron density states of the 
free electron gas is a direct measure of the energy level 
diagram of He I [l8| . The energy resolution of the spec- 
trum is 81meV (bin size of the histogram), sufficient to 
resolve the lowest lying triplet states in the range of 19.8- 
24.6eV as indicated. The structure that was observed in 
the continuum states of He I, as valence electrons asso- 
ciated to He + ranging between 24.6 and 40.8eV, is at- 
tributed to the molecular states of contaminating gases 
present in trace amounts. Note that we have an exact 
match of the electron energy distribution (a histogram 
of the raw data) with the known electron energy levels 
of He I, i.e. the observed spectrum is not corrected for 
an offset by the work function as in conventional surface 
Auger spectroscopy [9] . It is remarkable having observed 
the electronic eigen energies of helium in the simplified 
surface Auger experiment under evaluation. Apparently, 
the stationary fluctuations in the surface potential of a 
metal are directly connected to the eigen energies of the 
quantum system directly near the outer surface of the 
metal. In a first attempt to understand the thermody- 
namical aspects of this phenomenon, we have studied the 
effect under consideration in more detail using molecules 
(N2) instead. Unlike atoms, electrons in molecules can 
be treated as harmonic oscillators. Since the electronic 
energy levels in molecules are evenly and closely spaced 
(typically 200meV), the electron density of states can be 
approximated by a continuous distribution i.e. we ex- 
pect that the valence electrons of molecules behave as 
an idealized free electron gas in thermal equilibrium as 
suggested by eqn. [3J 



IV. OBSERVATIONS AND ANALYSIS 

In figure [3J we present a typical series of observed fluc- 
tuations in the electrode potential at fixed bias voltage 
when exposed to a nitrogen plasma. The stationary fluc- 
tuations that were studied were invoked only when the 
electrode was held in the region of the plasma jet were 
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(A) (B) (C) (D) 

FIG. 2: Schematic of the Auger de-excitation process of an exited state molecule near the surface of a negatively biased metal. 
Open circles represent hole states and black dots electron states. (A) An excited molecule, represented here as a two level 
(E a ,Eb) electron system, approaches a metal represented by its surface states. (B) At close range an electron, originally held at 
the surface, is transferred to the ground state of the molecule triggering the emission of the Auger electron into the continuum 
!• (C) The emission of multiple electrons at ambient conditions leads to the build up of space charge and elevates the Fermi 
level of the combined system. This process continuous until the Fermi level has increased by an amount of — qV — Et — E a . 
(D) If the excess of space charge is removed, the initial Fermi level of the system is maintained throughout the process. An 
event is then recorded as an increment in the voltage of the detector by an amount V. 



light emission occurs and completely disappeared when 
the electrode was held just outside this region. Next, the 
amplitude of the observed fluctuations increase with the 
electrode potential. The possible role of unintended elec- 
tronic interference could be ruled out by switching the gas 
discharge off while observing the afterglow of the nitro- 
gen plasma. The fluctuations could be observed during 
the decay of metastable nitrogen (see Appendix). There- 
fore, we can conclude that the detected fluctuations are 
no artifacts. 

The observed fluctuations at negative bias (figure [3p) 
relate to the emission of charge such that the voltage of 
the detector is increased. Under these conditions, elec- 
trons that were initially located at the surface of the 
electrode are transferred to the continuum by the Auger 
process (figure [5]). These events result in a temporary de- 
crease in the Fermi level and correspondingly to a tempo- 
rary decrease in the voltage as explained by the process 
discussed in figure [5] A similar process occurs at pos- 
itive bias. Electron energy distributions were obtained 
after building a histogram of the observed fluctuations. 
Typical results for zero, negative and positive bias are 
presented in figure O The energy resolution, defined by 
bin size of the histogram when sorting the recorded volt- 
ages, was 200meV. 

The autocorrelation of stationary fluctuations was 
found to decrease with increasing fluctuations in the volt- 
age. The correlation time of the signal when the plasma 
was turned 'off' (figure [HJ'V) exceeds I second. This was 
due to the finite resolution of the recording digital os- 



cilloscope. The data in figure [3] were acquired at a 
fixed setting of sensitivity of the oscilloscope. The de- 
tection resolution under these conditions was 20mV. The 
residual fluctuations in the voltage of figure [HJ'V (ther- 
mal noise) when measured at O.lmV resolution yields 
23±2mV (FWHM) (data not shown). Since the reso- 
lution of the measurements was nearly equal to the ther- 
mal noise level, the majority of successive readings of 
the voltage in figure [3]A were assigned the same value in 
the analog-to-digital conversion. This resulted in the ex- 
tended correlation time as observed in figure figure [SJ'V. 
Part of this effect, which is most pronounced at small 
amplitudes, is transferred to the fluctuations observed in 
figure [3)3. For signals with a large signal-to- noise ratio 
this effect is eliminated and an autocorrelation time is 
obtained that is typically less than the time between to 
successive readings. The correlation data in figures [3jU 
and dp satisfy K{t)/K(0) < 1/2tt for t > i.e. the 
readings in these series are uncorrelated. 



A. Coarse structure of the energy distribution 

We have recorded and analysed voltage fluctuations 
for a series of bias voltages in the range of -40 to +40V. 
The observed energy distributions of the two fractions 
were fitted to the energy distribution given by eqn. |3j to 
obtain the temperatures of both electron and hole en- 
sembles. At negative bias voltages (excess of electrons) 
we associate the energy distribution with a free electron 
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FIG. 3: Typical fluctuations (left) and autocorrelations (right) that were observed in the electrical potential of the Cu electrode 
exposed to a nitrogen plasma. (A) Base line which is recorded when the plasma source is turned off (B) Observed fluctuations 
when the plasma is turned on for the case of a zero biased detector, (C) a negative biased and (D) a positive biased detector. 
The autocorrelation time decreases with increasing fluctuations in the voltage. When the plasma source is activated the 
correlation time is reduced to less than 10 ms. For a biased detector the correlation time is typically less than the time 
between successive measurements. Under these conditions voltage readings are uncorrelated and readings can be considered as 
independent outcomes of a stochastic process. 



gas. This is illustrated in more detail in figure [5)3. Wc 
distinguish the elastic electron-electron scattering pro- 
cess (e-e) that occurs for Q < and E > from the 
inelastic hole-electron recombination process (h-e) that 
occurs for E < 0. Likewise, for a positive bias (excess of 
holes) we associate the energy distribution to a free hole 
gas (h-h) for Q > and E < and the decaying tail for 
E > that is the result of electron- hole (e-h) interactions 
(figure [5p). The energy distribution at zero bias is the 
composite of a residual free hole gas and a residual free 
electron gas. We suggest that the observed energy fluc- 
tuations at zero bias (Q = 0) is the direct result of two 
recombination processes: for E > we observe the en- 
ergy distribution of holes decaying in a free hole electron 
density (e-h) and for i?<0we observe the energy distri- 
bution of an electron ensemble decaying in an free hole 
density (h-e). In order to investigate these assumptions 
we have made an estimate of the lifetime of electrons 
and holes. We associate the energy distribution of an 
ensemble to the single particle lifetime using the Laplace 
transform. Under these assumptions the width of the ob- 
served energy distribution (AE) is a measure of the rate 



constant Q of the recombination process according to 

AE = fiF = - = kT. (4) 
r 

The recombination lifetime of electrons and holes un- 
der the distinct conditions discussed in this work are es- 
timated in table HI The lifetime of an electron at charge 
neutral conditions (zero bias) is 10.8±0.5 fs and lifetime 
of a hole 4.4±0.2 fs. The lifetime of a hole in an electron 
density yields 14±2 fs whereas the lifetime of an electron 
in an hole density yields 6.3±0.2 fs. 

The behavior of the temperature of the electron gas 
and the hole gas was analyzed in more detail (figure [6]). 
The temperatures of the hole gas and electron gas were 
determined (see Appendix) for two series at 15 different 
bias voltages ranging from -30 to +50V. The two data se- 
ries shown in figure [6] are the average of two independent 
measurements: the temperature at each bias voltage is 
a two point weighted average. The ensemble tempera- 
ture of the hole gas at negative bias is nearly a constant 
at the level A;T=47±6 meV. At positive bias the tem- 
perature inclines as a function of the bias voltage. The 
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FIG. 4: Observed photon- and electron energy spectrum in a helium gas discharge operated at 300K and 1 bar. (A) The 
primary content of the plasma is singly excited triplet helium as shown by the light emission spectrum. (B) The electron energy 
distribution that was obtained after analysis of the fluctuations in the static potential of the Cu electrode. As the valence 
electron in excited state helium can be considered as a free electron, the measured electron energy distribution coincides with 
the energy level diagram (density of states) of He I Note that no electron states were detected in the range of 0-19. 8eV. 
In the range of 19.8 to 24.6eV the lowest lying triplet states of He I were resolved and ascribed as indicated. The structure in 
the observed continuum states, ranging from the ionisation potential of He I at 24.6eV up to the threshold energy of He II at 
40.8eV, is attributed to molecular states of trace gases other than helium. 



temperature of the electron gas is nearly a constant, fixed 
to a level fcT=159±4 meV at negative bias. The expla- 
nation of latter behavior originates from the fact that 
the temperature of a degenerate gas is independent from 
the Fermi level. Actually, the chemical potential is inde- 



pendent of the temperature 0, HI ■ This is different for a 
gas at low densities (degenerate case). Upon raising the 
Fermi level under these conditions the kinetic energy of 
particles is continuously redistributed by elastic collisions 
within a range of A ~ kT, thus leading to an increase in 
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FIG. 5: Energy distributions derived from the observed fluctuations in the static potential of a nitrogen plasma. The energy 
distributions (data points) are fitted (solid lines) according to eqn. [3] (A) The energy distribution at zero bias is the resultant 
of two components which were separated by shifting the Fermi level of the system. (B) At negative bias electrons are added to 
the system. The energy distribution is observed is a free electron gas (e-e) for E > and a minor tail for E < which is due to 
the recombination process of an hole in an electron density (h-e). The inset at a log scale demonstrates the proper asymptotic 
behavior of the model. (C) At positive bias holes are added to the system and a free hole gas (h-h) is observed for E < and a 
minor tail for E > which is the recombination process of a electron in a hole density (e-h) . In this paper we suggest that the 
energy distribution that is observed at zero bias is the result of two processes: for E < it is an hole decaying in an electron 
density and for E > it is an electron decaying in a hole density. 



the temperature when particles are added. For a single 

component gas the temperature dependence of an degen- ^ 

erate gas was estimated by the elastic scattering rate of T = 1/r = -p av. (5) 

binary collisions Q: 

where p Q is the density of scatterers, a is the cross- 
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B. Molecular structure 

At high energy resolution, sharp, peaked structures 
were observed that occur at specific energy levels similar 
to the case of helium. The data shown in figure[7]is a typi- 
cal energy distribution of a hole gas for a nitrogen plasma 
detected at a resolution of 25 meV. The peaked structure 
in figure [7] resembles a molecular vibrational structure. If 
we account for the specific experimental conditions that 
were imposed, we deduce that these structures probably 
originate from redox reactions of exited state molecular 
nitrogen at the charged Cu electrode. This structure was 
also observed in the energy distribution of the electron 
gas. The observed free electron states coincide with the 
vibrational states of iV2(X 1 E+,^ = 0, ..,5) well below 
the threshold energy of iVrT at 1.6eV. 



FIG. 6: Temperatures of the observed free electron and free 
hole ensembles in a nitrogen plasma as a function of the bias 
voltage. The model for the h-h process is (solid line) is a 
straightforward extension of eqn. [5] and relates to an ideal 
free fermion gas at high temperatures. 



section of the process and v the relative velocity of the 
particles. The factor of | is to account for the sym- 
metry of the collision system (e-e or h-h). The relative 
velocity of two particles is given by v = y/2Ef/m. The 
density multiplied by the cross section in eqn. [5] can be 
expressed in terms of y/Ej (see eqn. [7]) . At high tem- 
perature (low density) conditions the width of the en- 
ergy distribution of identical particles is then given by 
HT = kT = 0.391E 1 /. Note, that this result is nearly 
identical to the thermal energy per particle given by 
kT = ~Ef in case of a free degenerate fermion gas Q. 
In our experiments the temperature of the free hole gas 
is assumed to present the equilibrium of a cold fraction 
(hole states in the ground state of N2) with a thermal 
energy fcT=25meV (T=300K) and a hot fraction (hole 
states in iV|) at thermal energy kT = 0.391E 1 /. In case 
of a diluted hot fraction, the equilibrium temperature of 
the mixture is a linear function in the number of parti- 
cles present in the hot fraction (ideal gas). In this work 
we have fitted the temperature dependence of the hole 
data in figure [5] to a linear function in the Fermi level 
to obtain a proportionality factor (experimental value) 
according to kT = (0.019 ± O.OOl)^/. 

Typical free electron and free hole densities that oc- 
cur in a nitrogen plasma were estimated using eqn. [5] 
at the corresponding temperatures listed in table |T] un- 
der conditions of zero bias. The electron density in 
the experiments evaluated for A:T e =149±8 meV yield 
p e = 2.6 ± 0.2 10 22 cm -3 and the hole density ph — 
0.68 ± 0.05 1 22 cm- 3 for kT h =6l±3 meV. 



V. DISCUSSION 

A. Electron emission process 

In presence of a metal surface the decay of excited 
state nitrogen is through the process of surface Auger 
effect. We have carefully checked whether the observed 
fluctuations as shown figure [3] were actually due to the 
activated nitrogen at the electrode under investigation 
as a result of the plasma jet. More specifically we have 
canceled the option that the emissions are the result of 
field emissions at high electric field strengths (corona 
discharges). In absence of the plasma, field emissions 
could only be observed after increasing the electric field 
strength at the electrode surface to a level in excess 
of 30 kV/cm. This was achieved after switching to a 
bias voltage of +(-)15kV using EMCO DX250(N) high 
voltage power supplies. The observed energy distribu- 
tion according to eqn. [3] is certainly different from elec- 
trons that are the result from high field field emissions 
d, |H [H, H3|. The energy distribution in field emis- 
sion follows the Fowler-Nordheim law Q as a result of 
quantum tunneling through the cathode-anode barrier. 
Secondly emissions that have been observed in our work 
occur in the energy range of to 20eV, which is one 
order of magnitude less than the threshold reported in 



B. Fluctuations 

Random fluctuations in energy and the relaxation of 
an excited state system to its ground state are subject 
to the fundamental thermodynamic aspects of finding 
the system in any of the accessible microscopic states 
[25| . Under stationary conditions, energy fluctuations 
obey the Wiener-Khinchine relations Q- The correla- 
tion K(s) of random temporal fluctuations in the en- 
ergy E(t) of an ensemble is given by the self convolu- 
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tianK(s) = (E(t)E(t + s)) 0. This correlation function 
is the Fourier transform of the spectral density of E{t) 
which is dn(E)/dE, if we assume that the electronic en- 
ergy of a quantum system satisfies E = fko. The exact 
relationship is given by 



dn(E) 
dE 



= F[K(s)] 



(6) 



The Wiener-Kinchine relations imply that the energy 
spectrum of an ensemble of electrons is directly obtained 
after taking the inverse Fourier transform of the auto- 
correlation of the time series qV(i). A constraint is that 
the voltage V(t) has to be recorded using a detector with 
a bandwidth exceeding qV/h. For the experiments un- 
der consideration this yields a typical bandwidth of 10 
PHz=10 16 Hz. Such devices do not exist. Therefore, a 
direct Fourier transformation using E(t) to obtain the 
density of states is only relevant from an academic point 
of view. However, employing a detector of such a band- 
width is not a requisite to determine the spectral density. 

A quantum system oscillates between its eigenstates 
(eqn. QJ at frequencies on the order of u>i = Ei/H until 
the event of a measurement (e.g. sampling the voltage 
of an electric system) when its eigenvalue Ei = —qVi 
is read. During the period that a reading is taken, 
the system under investigation is in a definite eigenstate 
which is maintained (in our case by a sample and hold 
circuit) until a successive reading is performed. Only 
after performing a successive measurement, the eigen- 
state will have changed randomly. This random change 
can be understood by analysing the time dependent re- 
sponse of a detector for signals that change at frequencies 
larger than the bandwidth. The time dependent volt- 
age can be decomposed into its harmonics according to 
V(t) = Visin(uit) + V 2 sin{uj 2 t) + .... Samples of V(t) 
that are recorded at a rate ui s << u>i yields the collec- 
tion of readings of V(t) at random phases. Under these 
conditions the likelyhood of an observation at level V is 
given by the probability density function of V(t) yielding 
D{V) = l/(n,/V? - V 2 ). Note that D{V) is just a lim- 
iting case of the Dirac delta distribution 5(Vi — V). Thus, 
the outcome of a series of measurements are static, sin- 
gle valued (quantised) readings. The expectation value 
for the energy of a free electron in an eigenstate is given 
by (Ei(t)) = (Vi(0)|ff|ViW) = (E iS in(Eit/h)). Using 
Ei = qVi we obtain u>i = Ei/% = qVi/h. Thus, a free 
electron with energy Ei gives rise to a harmonic signal 
of frequency a;, = qVi/h and amplitude Vi — Ei/q. As in 
the present work oj, is much larger than the bandwidth 
of the detector, the measurement yields the recording 
of a static value: (Ei(t)) = 8{Ei — E). In conclusion: 
the detection of the eigen energies of electrons yields the 
random observation of the amplitudes of high frequent 
harmonics, manifested as fluctuations in the static po- 
tential. 

This means that a spectrum of electronic eigen en- 
ergies can be obtained by the method employed in the 



manuscript, consistent with eqn. [6] The only require- 
ment is that the response time of the voltage probe is 
short compared to the time between successive measure- 
ments. In the experiments that have been presented 
this constraint was met since the response time of our 
probe which is RC=lpF-8.2Mf7=8.2/is (figure [TJ, the 
time needed to charge the probe, is much less than the 
time (=1 ms) between successive measurements. 

The probability that a transition occurs from a state 
with its energy in the interval {E, E + dE) to the ground 
state, is proportional to the density of states at energy E 
by virtue of Fermi's Golden rule [JOjIHl- It is therefore 
sufficient to determine the relative frequencies of uncor- 
rected readings of E(t) in the interval (E, E + dE) for a 
large number of measurements. The energy distribution 
dn(E)/dE is readily obtained by building a histogram 
of a sequence of readings of E(t). When voltage read- 
ings are uncorrelated (figure [3]) the density of states of a 
system is obtained by sampling and counting random val- 
ues of the energy as explored in this paper. Note that a 
complete reconstruction of the electronic wave function, 
a procedure of paramount interest to quantum comput- 
ing [l9l |. cannot be performed exclusively on the basis of 
known energy eigenstates. A complete reconstruction of 
the wave function requires, in addition, the determina- 
tion of the relative phases between the different eigen- 
states. These phases could not be determined as a result 
of the limited sample rate (bandwidth) of the detection 
system in present work. 

The peaked structures that were observed in all of our 
experiments was analysed in more detail using the typ- 
ical distribution of hole energies shown in figure [Jj As 
it was confirmed by lifetime measurements and photo- 
spectroscopy that the nitrogen jet primarily consist of 
excited state molecular nitrogen, the electron spectrum 
in figure [7j should contain part of the Auger spectrum of 
exited state molecular nitrogen at the Cu electrode sur- 
face @ . In figure [7J we have marked the onset energies 
(N^f = 0)) of characteristic vibrational states of the dif- 
ferent molecular symmetries indicated . We conclude 
that the energy distribution of the observed free hole gas 
is partly associated to the reduction states of molecular 
nitrogen. 



C. Interaction of electrons and holes 

Our aim in present work was the decomposition of a 
single redox reaction into its two fundamental half reac- 
tions. One reaction involving the addition of electrons 
(reduction of the free hole gas) and a second reaction 
involving the addition of holes (oxidation of the free elec- 
tron gas). It is well known that an oxidation and its 
corresponding reduction process at thermal equilibrium 
can proceed at different energies. An example is photo- 
synthesis by micro-organisms (26j where the energy dif- 
ference between two different redox states is balanced by 
the photon energy. The fundamental interaction process 
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FIG. 7: High resolution plot of the energy distribution of a free hole gas observed in a nitrogen plasma, revealing residual 
structure. The peaked structure coincides with known reduction states of excited state nitrogen that originates from the 
surface Auger process at the metal target. The molecular symmetry and electronic energies of the lowest lying N^fv = 0) 
states have been indicated. The ground state level of the molecular ion N+ lies at 15.6eV. 



of an electron and a hole is recombination [8J, |27j. In 
the case of a vacuum, recombination is proceeds by an- 
nihilation where an electron-hole pair recombines under 
emission of two 511keV photons (conservation of energy). 
Under the experimental conditions in this work, recom- 
bination of free electrons and holes that are associated to 
exited state metastable nitrogen occurs through sponta- 
neous decay by a singlet-triplet transition. The (kinetic) 
energy and angular momentum in the process of spon- 
taneous decay is conserved by emission of 3eV photons 
(see Appendix) whereas the rest mass is balanced by the 
molecule, acting as a third body itself. In the presence 
of a surface, recombination may occur through the sur- 
face Auger process were the vibrational excitation is lost 
by interaction with the electrode (figure [2]). The predic- 
tion of the lifetime of a hole (electron) in an electron(hole) 
density p for the latter process is estimated using binary 
collision theory Q according to 



T = l/r = pocrv (7) 

with a — nr 2 = 7r(47rp /3)~ 2 / 3 the cross-section of 
the scattering process, v the relative velocity of the elec- 
tron with respect to the hole given by v 2 — SkT/m 
(m=4.5 10 -31 kg is the reduced mass of the electron- 
hole system) and p the density of scatterers. The den- 
sity of electrons and holes is estimated by the free elec- 
tron (hole) density [f| which yields p = Zp m /Mm p with 
Z=3 (thus accounting for a trivalent N — N bond), M=28 
and m p = 1.68 10 _27 kg and p m = 1.25 kgm~ 3 for a ni- 
trogen environment at 1 bar and 300K. The calculated 
lifetimes according to eqn. evaluated at the temper- 
atures in the experiments, are given in table |TJ These 
theoretical predictions for the recombination lifetime of 
holes and electrons are in fair agreement with the experi- 
mental values if we take into account the crude estimates 
that have been made in adapting a free electron model. 
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TABLE I: Recombination lifetimes of free electrons and free 
holes in a nitrogen plasma under the experimental conditions 
that were distinghuised in this work. Experimental values for 
the lifetimes are according to eqn. U using the experimental 
ensemble temperatures. Theoretical estimates of the lifetimes 
are according to eqn. [7] Experimental values for the lifetimes 
are according to eqn. fusing the experimental ensemble tem- 
peratures. Theoretical estimates of the lifetimes are according 
to eqn. [7] 

Experimental Ensemble temp. Lifetime (fs) 



conditions kT (meV) experiment calculation 



negative bias 


47 ± 6 


14 ± 2 (h-e) 


9 


zero bias 


61 ± 3 


10.8 ± 0.5 (h-e) 


7.5 




149 ± 8 


4.4 ± 0.2 (e-h) 


4.8 


positive bias 


105 ± 4 


6.3 ± 0.2 (e-h) 


5.8 



We conclude that the interactions of the observed free 
electron and free hole gases may proceed through recom- 
bination and can be understood quantitatively. 

Based on the previous analysis we explore the sugges- 
tion that electrons and holes are different entities i.e. at 
thermal equilibrium hole states are not equivalent to elec- 
tron vacancies. Dirac's equations for the free electron 
incorporates the electromagnetic field (photons) and has 
four solutions. Except for two different spin states two 
types of solutions are known. One of the solution repre- 
sents the free electron in presence of the electromagnetic 
field. The second solution represents the conjugate state 
which has no physical reality. However, it is generally rec- 
ognized that this conjugate state (with negative energy 
solutions) are of significance both physically and mathe- 
matically [27|. This fundamental aspect on the behavior 
of electron and hole ensembles is treated by the theory of 
Quantum Electro Dynamics (QED). Despite this fact, it 
is highly controversial [28| to assert that free hole gases 
may represent this conjugate state since it is commonly 
assumed that a direct observation of these states is not 
possible [13, HI]. The conclusion that is justified, is the 
distinction between the energy distributions of free parti- 
cles with negative charge and positive charge. Hence, we 
have to assign different lifetimes to the processes of elas- 
tic scaterring and recombination of electrons and holes. 
This anomalous behavior of particles of opposite charge is 
remarkable since replacing q by -q (charge conjugation) 
in the experiment under consideration should not have 
affected the scattering rates as observed. 

We observed that the energy distribution of the 
observed free hole gas may be partly associated to the 
reduction states of molecular nitrogen. However, the 
energy distribution can not exclusively originate from 



the surrounding nitrogen. If this were the case, we would 
have detected a pure molecular state. It is well known 
that the wave function and eigen energy in eqns.[T]and[2] 
account for the combined system of detector+molecule 
i.e. it is the eigen energy of a mixed state [nj. A detailed 
study of this aspect requires careful consideration of 
the inherent invasive nature of quantum detection [lfj. 
This topic is of considerable interest in ongoing research 
fl9L |29| but lies beyond the scope of this paper. 



VI. CONCLUSIONS 

We have presented an in situ method to detect the elec- 
tronic eigen energies in the afterglow of gas discharges at 
ambient conditions. The principle of operation relies on 
the observation and analysis of stochastic fluctuations in 
the static potential of a single electrode. By using this 
method we could derive the energy distribution of elec- 
trons and holes present near the surface of the electrode. 
In a helium discharge, the electron spectrum coincides 
with the electronic states in He I. A free hole gas was 
observed in a environment of metastable nitrogen. The 
interaction with the electrode and excited state molecules 
surface was through the surface Auger process. The ex- 
change of charge is directed through omnipresent free 
electron and free hole gases. Recombination of electrons 
and holes was found to be the fundamental underlying 
process of reduction and oxidation. At ambient condi- 
tions the lifetime of a hole in an electron density yields 
4.4±0.2 fs whereas the lifetime of an electron in an hole 
density yields 10.8±0.5 fs. These measured lifetimes are 
in fair agreement with theoretical estimates of the pro- 
cesses which yield respectively 4.8 and 7.5 fs. Based on 
the observed asymmetry in the energy distributions of 
electron and hole ensembles we suggest that under con- 
ditions of thermal equilibrium the behavior of a hole is 
not equivalent to that of an electron vacancy. 
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Appendix 
A. Metastable molecular nitrogen 

A plasma generated in a nitrogen environment at am- 
bient conditions using a high voltage discharge contains 
a considerable amount of metastable nitrogen. Molec- 
ular metastable nitrogen A^(A 3 £+), with an internal 
energy of 6.2eV and an estimated radiative lifetime of 
1.9 sec [l6j | , w as previously observed at ambient con- 
ditions by PH. We have verified that N 2 {A 3 T,+) was 
produced by our specific source by analysing the after- 
glow of the plasma using emission spectroscopy and life- 
time measurements. Emission spectra were recorded per- 
pendicular to the cold plasma jet using a fiber probe 
located 20mm above the nozzle. Light was analysed 
using an Avaspec 2048-2 UV-VIS spectrometer system 
supplied with a FC-UV400-2 fiber, UA-178-1100 grat- 
ing and a Sony 2048 DUV detector board. We observed 
light emissions from the A 3 Y<+(v = 1, 2,3.1^ 1S g(^ = 
10, 11, 12, 13) system (Vegard-Kaplan) [lj, |l6|. In addi- 
tion we verified the lifetime of the excited state molecules 
present in the jet by Residence Time Distribution anal- 
ysis. After connecting a drift tube (Tricolair hose 
10x16mm, 3m in length) to the nozzle we found that 
the light emission of the afterglow was extended over the 
entire length of the tube. The lifetime of the metastable 
fraction was analysed based on the survival of fraction 
exited state species detected at the end of the drift tube 
as a function of the Residence Time (RTD analysis). The 
excited state molecules were produced in bunches by pe- 
riodically switching the high voltage of the discharge off- 
on-off with a 10% duty cycle. We detected the relative 
numbers and the distribution of arrival times of the sur- 
viving excited state fraction using a Faraday cup (alu- 
minia grid 20x20mm, 3mm mesh, 70% open) located at 
the outlet of the tube. The emission current was detected 
using a Keithley 610C electrometer as a function of the 
time delay between the pulsing of the discharge and the 
arrival of excited state molecules at the detector. The gas 
flow was adjusted in the range of 0.1-0.3 L/s at laminair 
air flow conditions to control the residence time of excited 
state nitrogen in the drift tube. Typical residence time 
distributions and a corresponding emission spectrum are 
plotted in figure^ The spectroscopic data in conjunction 
with the observed lifetime of the excited state nitrogen 
is attributed to the presence of iV 2 (A 3 E+). 

Cold plasma's based on electric discharges in nitro- 
gen (figures [1] and [H]) that operate under ambient con- 
ditions have been studied since 1916 fl2l. [l3l [lil 
Spectroscopic studies in the afterglow [16j | of nitrogen 
have revealed the presence of metastable molecular ni- 
trogen. We could confirm that the giant afterglow [l3[ 
in nitrogen is due to metastable N 2 by combining emis- 
sion spectroscopy with lifetime measurements (figure [9| . 
It appears that A^A), which has an internal energy of 
-EU=6.2eV, can exist at concentration of approximately 
1 ppm in a thermal equilibrium with ground state N 2 
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FIG. 8: The jet of excited state nitrogen (cold plasma) is 
based on a high voltage discharge operated at ambient con- 
ditions according to 

0, El El- The throughput of the car- 
rier gas is typically 0.2 L/s. The purity of the nitrogen was 
99.5±0.2% which was determined using a Varian Micro MC- 
CP-2002 gaschromatography. The temperature of the cold 
plasma flame, measured using a mercury thermometer, was 
40±2°C. 



at ambient temperatures (figures [8] and [9]) . This is re- 
markable, as under these experimental conditions A^(A) 
is able to survive more than 10 9 exchange reactions 
N 2 (A) + N 2 -* N 2 + N 2 (A) before detection. One could 
argue, in order for the second law of thermodynamics 
to hold, that only a fraction of magnitude riA/n — 
e -EjE A = e -M7-E.4 = 2.0TCr 100 should be found in the 
present experiments. However, this argument is incor- 
rect. The reason why an extreme afterglow can be ob- 
served, even at ambient conditions, is two fold. First, the 
lowest exited state A 3 Y;+ is a triplet state whereas the 
ground state is a singlet. Radiative transitions between 
triplet and singlet are forbidden by spin selection rules 
and therefore unlikely to occur. Secondly, N 2 (A) is sta- 
ble in a N 2 environment as a result of the symmetry of 
the exchange reaction N 2 (A) + N 2 — > N 2 + N 2 (A) which 
additionally excludes de-excitation by means of two body 
interactions. Similar to the nuclear Overhauser effect Q 
observed in spin polarized gases [13, Hl[ the relaxation 
of N 2 (A) is only by virtue of higher order de-excitation 
processes: the system is inherently non-stationary. Thus, 
despite that the fraction of N 2 (A) has completely ther- 
malised, no true stationary equilibrium is reached in the 
afterglow. Thermodynamical laws simply do not hold 
[25l ] under these 'extreme' conditions, thereby resolving 
the paradox. We conclude that N 2 (A) is stable in a N 2 
environment and that its decay is primarily through ra- 
diative emissions of forbidden light (Vegard-Kaplan) [lj| . 



B. Voltage probe 

In order to reduce the load capacitance and simultane- 
ously increase the input impedance of our voltage probe 
(figure Q} we used a SMD resistor (MMB0207 Vishay 



BC components) of 8.2MOhm//lpF (the residual stray 
capacitance across the SMD resistor yields ~lpF) in sc- 
ries with the 1 MOhm//12pF input impedance of the 
recording oscilloscope (Tektronix TDS-3052B). Note that 
this probe operates according to the scheme of a nearly 
matched, high frequency voltage probe [TtJ- However, 
a connection cable between the high impedance resis- 
tor and the IMOhm input resistance was avoided in 
our set-up. In this way the load capacity could be re- 
duced to a pF level. This high impedance, low capacity 
probe allows the detection of voltage fluctuations at rel- 
atively small charging currents. The voltage readings of 
our probe were calibrated against a 30.0±0.1 VDC refer- 
ence. Alternatively, we used a Tektronix P6139A (10:1) 
voltage probe. The P6139A has an input impedance of 
10MOhm//8.0pF. Its relative high load capacitance is 
due to the compensation of the 0.7m long cable between 
the tip and the connector. The free hole and free electron 
gases reported in this work were observed using cither 
the 8.2MOhm//lpF probe or the P6139A probe. The 
voltage reference of the recording oscilloscope and probe 
was biased using an external DC power supply (Hewlett 
Packard 412B) as indicated in figured] 

C. Data reduction and fitting routine 

The fluctuations at a fixed bias level were analysed 
by building a histogram (bin sizes 200 meV or 25 meV) 
of the observed voltages. It should be obvious that the 
choice of the bin size when building a histogram affects 
the resolution of the observed energy distribution. A 
time series V(t) typically contains 10.000 readings and 
series were recorded at a rate of 1kHz. Each data point 
V(t) was converted to a sample of the energy Ei using 
Ei — —qV(ti). The energy distributions thus obtained 
were fitted by means of a non-linear least-squares routine 
using the function: 

ME) = {B e VE [exp [p e (E - E s )\ + If 1 E > 
dE \B h y/^Eexp[(3 h E] E<0 

where (3 e = fcT e , E s and an overall scaling factor B e are 
the fitting parameters for the electron gas and /3h = kT^ 
and an overall scaling factor arc the fitting parameters 
for the hole gas. To account for a mutual offset in the 
energy scale, the fit function in eqn. [5] was allowed to 
shift (using an additional fitting parameter). Typically, 
the offset compensation yields approximately 120 meV. 
By using this fitting routine estimates were obtained for 
the temperature of the electron gas, the temperature of 
the hole gas and for the Fermi energy of the electron gas. 
Based on the information contained in all our data it 
is concluded that the hole gas corresponds to a fermion 
gas in the limiting case of high temperatures and low 
density. The Fermi-Dirac distribution factor in the fitting 
routine was therefore reduced to a Maxwell-Boltzmann 
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FIG. 9: The molecular content of the nitrogen plasma jet was analysed by emission spectroscopy and lifetime measurements. 
(A) The characteristic forbidden A 3 £+0 = 1, 2, 3)-X 1 E+(f = 10,11,12,13) Vegard-Kaplan emissions that were observed in 
the afterglow. (B) The lifetime of the excited state molecules, that was studied by Residence Time Distribution analysis using 
a 3 meter long drift tube, yields a lifetime on the order of 1 second. 



distribution. This behavior is different for the electron 
gas. The electron energy distribution could be fitted only 
when fully accounting for the Fermi-Dirac statistics: the 
electron energy distribution has a non-zero shift in the 
Fermi energy. At negative bias E s typically yields 400 
meV at a temperature of fcT e =150 meV, implying that 



we deal with the non-degenerate case (high density low 
temperature limit) of a free electron gas. The estimates 
of the fitting parameters obtained by this fitting routine 
reproduces well if data are acquired for a fixed position 
of the electrode (25±5 mm) in the plasma flame. 



